Dynamical r-process studies within the neutrino-driven wind scenario and its 

sensitivity to the nuclear physics input 
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We use results from long-time core-collapse supernovae simulations to investigate the impact of 
the late time evolution of the ejecta and of the nuclear physics input on the calculated r-process 
abundances. Based on the latest hydrodynamical simulations, heavy r-process elements cannot be 
synthesized in the neutrino-driven winds that follow the supernova explosion. However, by artificially 
increasing the wind entropy, elements up to A — 195 can be made. In this way one can reproduce 
the typical behavior of high-entropy ejecta where the r-process is expected to occur. We identify 
which nuclear physics input is more important depending on the dynamical evolution of the ejecta. 
When the evolution proceeds at high temperatures (hot r-process), an (n, 7) «^ (7,1) equilibrium is 
reached. While at low temperature (cold r-process) there is a competition between neutron captures 
and beta decays. In the first phase of the r-process, while enough neutrons are available, the most 
relevant nuclear physics input are the nuclear masses for the hot r-process and the neutron capture 
and beta-decay rates for the cold r-process. At the end of this phase, the abundances follow a steady 
beta flow for the hot r-process and a steady flow of neutron captures and beta decays for the cold 
r-process. After neutrons are almost exhausted, matter decays to stability and our results show 
that in both cases neutron captures are key for determining the final abundances, the position of 
the r-process peaks, and the formation of the rare-earth peak. In all the cases studied, we find that 
the freeze out occurs in a timescale of several seconds. 

PACS numbers: 26.30.-k, 26.30.Hj, 26.50.+X, 97.60.Bw 



I. INTRODUCTION 



The synthesis of heavy elements by the rapid neutron 
capture process (r-process) is a fascinating, long-standing 
problem which involves challenges in nuclear physics (ex- 
periment and theory), astrophysical simulations of explo- 
sive environments, and observations of metal-poor stellar 
atmospheres. The astrophysical site where the r-process 
produces half of the heavy elements has not yet been 
identified (see Ref. [l[ for a recent review). However, 
even if the astrophysical scenario were discovered, one 
still has to deal with nuclei far from stability for which 
no experimental data are available and theoretical pre- 
dictions are quite uncertain [2|. 

Galactic chemical evolution models 3, 4] indicate that 
heavy elements are most likely produced in core-collapse 
supernova outflows. After a core-collapse supernova 
explosion, a proto-neutron star forms and a baryonic 
neutrino-driven wind develops The matter expands 
at high velocity, which can become supersonic, and even- 
tually collides with the slow, early supernova ejecta re- 
sulting in a wind termination shock or reverse shock 
lioj . There are several nucleosynthesis process that occur 
or might occur in this environment: a-process [TlL fl2j . 
i^p-process fl3l - [l5j , and r-process [H, E3 ■ The produc- 
tion of heavy r-process elements (A > 130), requires a 
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high neutron-to-seed ratio. This can be achieved by the 
following conditions [l8l - [2l| : high entropy, fast expan- 
sions, or low electron fraction. However, these conditions 
are not yet realized in hydrodynamical simulations that 
follow the outflow evolution duringthe first seconds of 
the wind phase after the explosion [H, HoL |2^|. 

This manuscript aims to explore the sensitivity of the 
calculated r-process abundances to the combined effects 
of the long-time dynamical evolution and nuclear physics 
input providing a link between the behavior of nuclear 
masses far from stability and features in the final abun- 
dances. The impact of the nuclear physics input on the 
production of heavy elements has been explored in clas- 
sical r-process calculations in Ref. [23], [24| and with more 
dynamical but still parametric calculations (e.g., [25l — 
[23I )• However, most of the r-process studies [l| aimed to 
find the astrophysical conditions (density, temperature, 
electron fraction) that reproduce the solar abundances 
for a given nuclear physics input. There have been only 
few works ex plor ing the impact of different mass mod- 
els (e.g., [28l. l29l|). the effect of neutron captures when 
matter decays to stability [3014331 ] , and the importance of 
beta decays [H], HU . Our results contribute to improve 
the present understanding of how nuclear masses and 
neutron-capture cross sections determine the r-process 
abundances. The influence of the beta-decay rates will 
be studied in future work, here we only explore the effect 
of beta-delayed neutron emission. 

In this paper we use hydrodynamical trajectories from 
the neutrino-driven wind simulations of Ref. [|. As the 
conditions found in these trajectories do not allow for the 
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synthesis of heavy r-process elements [36| , we artificially 
increase the entropy by a factor two in order to produce 
the third r-process peak. This allows us to study the nu- 
cleosynthesis of heavy elements in a typical high-entropy 
neutrino-driven wind in a more consistent way than with 
fully parametric expansions [29| or with steady-state 
wind models (e.g. [2(1 HH), which cannot consistently 
explore the interaction of the wind with the slow super- 
nova ejecta. The possible influence of the wind termina- 
tion shock on the nucleosynthesis was already suggested 
in Refs. @, EH an< ^ analyzed in different works [37l - l4lj |. 
For the first time, we use here trajectories obtained in 
hydrodynamical simulations to perform a detailed inves- 
tigation of the effect on the r-process of the long-time 
dynamical evolution. We choose two different astrophys- 
ical evolutions, that cover the broad range of conditions 
found in the hydrodynamical simulations. Four different 
mass models are used to study the impact of the nuclear 
physics input on the calculated abundances. 

Our astrophysical and nuclear physics inputs are in- 
troduced in Sect. [II] We study the impact on the final 
abundances of the reverse shock and long-time dynamical 
evolution fSect. lHTA")) . The sensitivity of the abundances 
and of the r-process evolution to the mass model is ex- 
plored in Sect. IIIIBl In Sect. IIII CI we discuss the evo- 
lution of the abundances after the r-process freeze-out, 
that we define as the moment when the neutron-to-seed 
ratio drops below one. Our conclusions are summarized 
in Sect. E3 



II. METHODS 
A. Neutrino-driven wind and termination shock 

Our nucleosynthesis studies are based on a trajectory 
ejected at 8 s after bounce in an explosion of a 15 M Q 
progenitor (model M15-12-rl in [9j). This trajectory rep- 
resents a typical neutrino-driven wind, whose wind phase 
can be described by steady-state models [H, HH, IHJ . 
However, hydrodynamical simulations are required to 
study the long-time evolution when the supersonic wind 
collides with the slow-moving supernova ejecta resulting 
in a wind termination shock. 

Detailed description of the hydrodynamical simula- 
tions can be found in Refs. M, |43j. In these simulations, 
Newtonian hydrodynamics [43l . l44j with general relativis- 
tic corrections for the gravitational potential 45[ is com- 
bined with a simplified neutrino transport approxima- 
tion assuming Fermi-Dirac neutrino spectra 43[. This is 
computationally very efficient and reproduces the results 
of Boltzmann transport simulations qualitatively. The 
equation of state used in the simulations includes neu- 
trons, protons, alpha particles, and a representative nu- 
cleus ( 54 Mn) treated as non-relativistic Boltzmann gases 
in nuclear statistical equilibrium Q. The central part 
(p > 10 13 g/cm 3 ) of the proto-neutron star is removed 
from the computational domain and a Lagrangian inner 



boundary (placed below the neutrinosphere) describes 
the neutron star evolution. The neutrino cooling (i.e., 
neutrino energies and luminosities) and contraction of the 
proto-neutron star are parametrized at the inner bound- 
ary to account for possible uncertainties in the high den- 
sity equation of state. 

The neutron-to-seed ratio found in the simulations of 
Ref. [9( after freeze-out of charged-particle reactions is 
too low (Y n /Y scc a ~ 10 -2 ) to permit the formation of 
heavy r-process nuclei. Only elements with Z < 48 are 
produced, i.e. light element primary process (LEPP) ele- 
ments [36} . However, we can still use this trajectory for r- 
process studies, if the neutron-to-seed ratio is artificially 
increased by assuming a smaller initial electron fraction 
or by rising the entropy. The electron fraction is deter- 
mined by electron neutrino and antineutrino energies and 
luminosities (see e.g., Ref. [l8j]) and we keep it as given by 
the simulations (Y e w 0.47). Hereafter, all calculations 
are performed on the same trajectory with the density 
decreased by a factor of two [9(| overall to get also a fac- 
tor two higher entropies (S oc T 3 /p w 200 fc^/nuc) and 
thus high neutron-to-seed ratio (Y n /Y see d ~ 70). This 
is enough to produce nuclei around the A = 195 peak 
and mimics the hydrodynamical conditions of a neutrino- 
driven wind where the r-process does occur and of other 
astrophysical environments that involve ejection of high 
entropy matter. Therefore, it can be used as basis to 
study the combined influence on the abundances of the 
long-time dynamical evolution and of the nuclear physics 
input. 

Variations in the late evolution of the ejecta are ex- 
pected as shown in multidimensional simulations [46| . 
Although the neutrino-driven wind stays spherically sym- 
metric in absence of rotation, the interaction of the wind 
with the slow, early supernova ejecta and thus the result- 
ing reverse shock depends on the progenitor structure @ 
and on the anisotropic pressure distribution of the su- 
pernova ejecta, where the wind propagates through [46| . 
Since our trajectory corresponds to a spherically sym- 
metric simulation, we have modified it to account for 
the possible variation of the reverse shock radius and of 
the long-time evolution. The changes are done only af- 
ter charged-particle reactions freeze out to assure same 
initial conditions for the r-process. 

As we will show, the reverse shock has a non-negligible 
influence on nucleosynthesis. When the wind collides 
with the slow-moving ejecta, the expansion velocity drops 
as kinetic energy is transformed into internal energy, with 
the consequent increase in temperature and density. The 
density (p r s), temperature (T r s), and velocity (u rs ) of 
the shocked matter can be calculated with the Rankinc- 
Hugoniot conditions, corresponding to mass, momentum, 
and energy conservation: 
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here p, u, P, and e are the density, velocity, pressure, 
and specific internal energy, respectively. Quantities in 
the wind have the subscript "w" and in the shocked ma- 
terial, above the reverse shock, the subscript "rs". In 
addition to these equations, one needs an equation of 
state that relates pressure and energy. In the last equa- 
tion we take e ~ 3P/p which is a good approximations 
because the environment is radiation-dominated. The lhs 
of these equations is known from the wind, therefore we 
combine them to obtain two possible solutions for the 
matter velocity after the shock: 1) u IS = u w , no shock; 
2) w rs = u w /7 + 8/7P w /(p w )i w ). Once the velocity is 
known, the density and pressure of the shocked matter 
are calculated with Eqs. (fTaj). ([Tb|) . Any other thermo- 
dynamical variable, including temperature and entropy, 
can be obtained from an equation of state (EoS). Here, 
we use the Timmes EoS [47] . 

Once the conditions of the shocked matter are known, 
we need to describe their evolution. The density and ve- 
locity evolutions have to fulfill the condition of constant 
mass outflow (M = 4irr 2 vp). Two extreme expansions 
can be identified: 1) the velocity is constant, and the 
density thus decreases as r -2 ; 2) the density is constant 
and then it is the velocity that decreases as r~ 2 . The 
latter expansion was used in Ref. but it implies a 
decrease of the velocity down to a few m s _1 in about a 
second, while in full hydrodynamical simulations 0, [Io| , 
the velocities stay around 10 3 -10 4 km s _1 . Following 
these simulations, we prescribe an extrapolation for the 
evolution of matter which is between these two extreme 
cases. The density of the shocked matter stays constant 
and the velocity decreases as r~ 2 during one second. Af- 
terwards, we keep the velocity constant and the density 
decreases as r~ 2 . A similar extrapolation is obtained 
with the prescription used in Ref. [48| . 



B. Nucleosynthesis network and nuclear physics 
input 

We start our nucleosynthesis calculations at a temper- 
ature T = 10 GK where the composition is given by 
nuclear statistical equilibrium and is dominated by free 
neutrons and protons. The evolution of the composi- 
tion is followed with two reactions networks. During the 
seed formation we use an extended network that includes 
all possible charged-particle reactions. While for the r- 
process we take advantage of a faster network that only 
considers neutron capture, beta decay, photodissociation, 
alpha decay, and fission. However, fission reactions play 
a negligible role in the present calculations and will not 
be further discussed. 

The extended nuclear reaction network consists of 
3347 nuclei from neutron and protons to Europium. 
Neutral and charged-particle reactions are the same as 
in the REACLIB compilation used by 49]. For the 
weak- interaction rates (electron/positron captures and j3- 
decays) we use the rates of Refs. [5(1 [H| for nuclei with 



A < 45 and those of Refs. [5J, [53j for 45 < A < 65. 
Neutrino interactions are important during the seed for- 
mation since they control the amount of free neutrons 
and arc taken from Ref. (54[. When the temperature 
drops below T ~ 3 GK charg ed-p article reactions freeze 
out, i.e end of the a-process [ill. H^|. and the r-process 
phase, characterized by a domination of neutron cap- 
ture, begins. The subsequent evolution is followed by 
a r-process network of 5300 nuclei between Z = 14 and 
Z = 110. Since we are interested in the detailed evolution 
of matter when decays to stability, we had to improve 
the original network of D. Mocelj PhD (5||. This net- 
work was based on the algorithm suggested by Ref. [56[ 
which is also used in Refs. [26|, [2{|. Such algorithm as- 
sumes that the neutron abundance (Y n ) stays constant 
during a time step (see Appendix) . This is very efficient 
during the early r-process phase, but it becomes numer- 
ically unstable when the neutron-to-seed ratio drops be- 
low one and Y n decreases very fast. This problem can 
be cured using smaller time steps, which increases how- 
ever the computational time without completely remov- 
ing the numerical instabilities. In our updated network, 
the equation for the neutron abundance evolution (see 
Eq. IA.3[) is included and the resulting system of equa- 
tions is solved using a fully implicit scheme based on the 
Newton- Raphson method 1571 158| | and the sparse matrix 
solver package PARDISO [59f , 

Our nucleosynthesis calculations are based on four dif- 
ferent mass models and their consistently calculated neu- 
tron capture rates: the Finite Range Droplet Model 
(FRDM) [6(|, the quenched version of the Extended 
Thomas- Fermi with Strutinsky Integral (ETFSI-Q) [6lj| . 
the version 17 of the Hartree-Fock-Bogoliubov masses 
(HFB-17) [H, and the Dufio-Zuker mass formula 
For the FRDM and ETFSI-Q mass models the neutron- 
capture rates are taken from Ref. [641) t hat uses the 
statistical model code NON-SMOKER [6^. The HFB- 
17 neutron capture rates are taken from the Bruslib 
database 9l| and were computed with the statistical 
model code TALYS jgi]. For the Duflo-Zuker mass for- 
mula, we have evaluated the neutron-capture rates using 
the analytic approximation suggested in Ref. (6?J that 
reproduces the results of more sophisticated statistical 
model calculations [68, [6^ (see also Sect. MI Cl) . For the 
range of temperatures we are considering, the tempera- 
ture dependence of the (n, 7) rates can be neglected and 
we use thus values that correspond to a temperature of 
30 keV. The inverse (7, n) rates are obtained from de- 
tailed balance (see Eq. IA.2I) . 

We use theoretical beta-decay rates from Ref. [7(| 
supplemented by experimental data whenever available 
(NuDat 2 database [7l|). We realize that the beta decays 
should be calculated consistently with the mass model, 
however such calculations have not been performed so 
far. We plan to explore the sensitivity of r-process cal- 
culations to different beta-decay rates in future work. 
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III. RESULTS 



Wind termination shock and long-time 
evolution 



Here we analyze the impact of the reverse shock on 
the r-process abundances and dynamics. All calculations 
presented in this section are based on the ETFSI-Q mass 
model. 
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FIG. 1: (Color online) Temperature and density evolution of 
a mass element ejected at 8 s after bounce of the supernova 
explosion based on model M15-12-rl of Ref. [jj, see discussion 
in Sect. [TTJ The solid black line corresponds to the original 
trajectory from the supernova simulation labeled as "unmod- 
ified". Notice, that density have been divided by a factor 2 
to get higher neutron-to-seed ratio. The green dashed line 
represents a evolution with the reverse shock at temperature 
of 1 GK. For the evolution shown by the red dashed-dotted 
line the reverse shock was removed. 



Our nucleosynthesis study correspond to the trajecto- 
ries shown in Fig. [TJ The solid black line, labeled as 
"unmodified" , represents the trajectory from Ref. @ in- 
troduced in Sect. Ill Al The position of the reverse shock 
and the evolution after it are not modified, but the den- 
sity is overall reduced by a factor of two. In the dashed 



green line the reverse shock is assumed to be at tem- 
perature of 1 GK and the subsequent evolution is calcu- 
lated as described in Sect. Ill Al The dashed-dotted red 
line, labeled as "no rs" , reproduces a case without reverse 
shock, where matter expands without colliding with the 
slow, early supernova ejecta. In this case, we assume an 
adiabatic (constant entropy) expansion with constant ve- 
locity starting at the position of the reverse shock in the 
simulation. As the mass outflow is constant, the density 
decreases with radius as p oc r~ 2 . 
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FIG. 2: (Color online) Final abundances for the different 
evolutions of Fig. ft] compared to solar r-process abundances 
shown by dots. 

Figure [2] shows the final r-process abundances for 
the different trajectories together with solar r-process 
abundances (dots) [72j . None of the calculations re- 
produce the solar abundances around the second peak 
(A ~ 130), since we have chosen the conditions which 
produce mainly the third r-process peak (^4 ~ 195). A 
pattern similar to solar is expected from a superposition 
of different trajectories [26|, EH • Some of the deficiencies 
seen in Fig. [2] are due to the mass model and will be dis- 
cussed in the next section. However, there are features 
that depend on the dynamical evolution. In order to un- 
derstand the abundances under different late time evolu- 
tions (Fig.QJ, we look at the characteristic timescales for 
the r-process: neutron capture, photodissociation, and 
beta decay that are defined, respectively, as: 
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where N n is the neutron number density, (ctv}(z,a) 
neutron capture or (n,j) rate, A 7 (Z, A) the photodisso- 
ciations or (7, n) rate, and \p{Z, A) the beta-decay rate 
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The evolution of these timescales is shown in Fig. [3] for 
the trajectories labeled as "T rs = 1 GK" (left panel) and 
"unmodified" (right panel). The trajectory labeled "no 
rs" follows a behavior very similar to the unmodified one. 

For the trajectory with the reverse shock at T = 1 GK 
(Fig. [31 left panel), there is a competition between neu- 
tron captures and photodissociation that last until neu- 
trons are exhausted at around 1 s. During this evolu- 
tion under (n, 7) ^ (7, n) equilibrium, the beta-decay 
timescale is longer than the other two, i.e. T( ni yj = 
T (7,n) ^ r /3- I n the following, this kind of evolution is 
called hot r-process. When the reverse shock is at lower 
temperatures (Fig.O right panel), the photodissociation 
timescale becomes longer than the other two timescales 
once the temperature drops below ~ 0.5 GK. The subse- 
quent evolution proceeds by a competition between beta 
decay and neutron capture, i.e. T(„ j7 ) « Tg <C i~( 7i „). This 
evolution was already studied in Ref. [73[ and has been 
recently named as cold r-process and rn-process [4l| . 
We call this evolution cold r-process. 

The relevant nuclear physics input depends on the dy- 
namical evolution. The hot r-process proceeds initially 
in (n, 7) ^ (7, n) equilibrium, therefore the neutron sep- 
aration energy is the key quantity that determines the lo- 
cation of the r-process path, see Eq. (|A.4p . For the cold 
r-process the evolution proceeds by a competition be- 
tween neutron captures and beta decays. Therefore, the 
most relevant nuclear physics inputs are beta-decay and 
neutron-capture rates. Nuclear masses are also impor- 
tant because they enter in the calculation of both. The 
evolution after freeze out is dominated by beta decays 
and neutron captures for both cold and hot r-process. 

The final abundances shown in Fig. [^indicate that cold 
r-process calculations lead to broader peaks because the 
r-process proceeds farther away from stability. The hot 
r-process, which evolves in (71,7) (7,1.) equilibrium, 
results in a huge trough in the abundances before the 
third r-process peak (green line in Fig. [2]). This is due 
to the behavior of the neutron separation energy just 
before the N — 126 shell closure (see Fig. [5]) and will be 
discussed in the next section. 

After we have introduced the two possibilities for the 
r-process: hot and cold, the evolution of relevant quanti- 
ties (neutron density, neutron-to-seed ratio, and average 
neutron separation energy) will be explained. The av- 
erage neutron separation energy (shown in the bottom 
panel of Fig. 0} is defined as: 
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We consider first the hot r-process which is labeled as 
"T rs = 1 GK" in Figs. Q] [2J and U The average neu- 
tron separation energy (Fig. HI bottom panel) decreases 
very fast initially, as matter moves away from stability 
after charged-particle reactions freeze out. (S n ) comes 
to a minimum after 30 ms, because the matter flow has 
reached N — 82 shell closure. Here the abrupt drop of 
individual neutron separation energies and the high pho- 
todissociation rates prevent matter to move farther away 
from stability. Therefore, the r-process path moves to 
higher Z by beta decays and successive neutron captures, 
while the neutron number stays constant at N = 82. 
This is shown in Fig. [6] by the dots that mark the r- 
process path. Once the neutron separation energy for 
nuclei beyond A" = 82 becomes large enough, the flow 
of matter can continue moving towards heavier nuclei. 
The r-process path reaches the N — 126 shell closure at 
around 500 ms as indicated by the second minimum in 
(S n ). When matter starts to pass through the N = 126 
shell closure, neutrons are exhausted in our calculations 
and the matter decays to stability. 

The evolution of the average neutron separation energy 
for the other two trajectories (cold r-process) is very dif- 
ferent because photodissociation becomes negligible once 
the temperature is < 0.5 GK. The evolution proceeds 
by a competition of neutron captures and beta decays 
and this allows matter to move farther from stability 
compared to the hot r-process. Therefore, the average 
beta-decay lifetime becomes shorter which speeds up the 
flow of matter towards heavier nuclei and broadens the 
minimum in (S n ). The faster evolution leads to a more 
rapid decrease of the neutron-to-seed ratio (Fig. SI mid- 
dle panel) and therefore the r-process ends earlier than 
in the hot r-process. 

The case without reverse shock (red line in Fig. [1]) is 
extreme because the neutron density decreases initially 
very fast (Fig. |H upper panel). This leads to a drop of 
the neutron captures, which explains the high values of 
the neutron density and neutron-to-sccd ratio that are 
maintained at later times. Having a large neutron-to- 
seed ratio at late times allows a continuation of neutron 
captures, even after several seconds. Therefore, the peak 
at A = 195 is shifted towards higher mass numbers as 
shown in Fig. [5] (see also discussion in Ref. |3l|). In the 
hot r-process, neutron captures after freeze-out also lead 
to the shift of the third peak, even when the r-process 
path at freeze-out is located at a neutron separation en- 
ergy of S n = 2.8 MeV, similar to the one used in the clas- 
sical r-process calculations of Ref. (23|, where the peak 
in the final abundances is obtained at the right position. 



where S n (Z, A) is the neutron separation energy of a nu- 
cleus with mass number A and charge number Z, and 
Y(Z,A) is its abundance. The neutron separation en- 
ergy is smaller for nuclei far from stability, since their 
neutrons become less bound. Therefore, we can use this 
quantity to study the evolution of the r-process. 



B. Influence of the mass model 

Nuclear masses are a key nuclear physics input for r- 
process calculations as they determine the energy thresh- 
olds for all relevant reactions: neutron capture, photodis- 
sociation, and beta decay. In this section we present 
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FIG. 3: Evolution of relevant timescales for the trajectories shown in Fig. [1] and labeled as "T rs = 1 GK" (left panel) and 
"unmodified" (right panel). Solid, dashed, and dotted lines represent the neutron capture, beta decay, and photodissociation 
timescales, respectively. 



results based on the four mass models introduced in 
Sect. Ill Bl and link features in the abundances with the 
behaviour of the two neutron separation energy (S^n). In 
the r-process evolution, two phases can be distinguished 
depending on whether the neutron-to-seed ratio is larger 
or smaller than one, i.e. before and after freeze-out. The 
abundances at freeze-out (shown in Fig. [5] by black lines) 
contain the information of the pre-freeze-out phase when 
there are enough neutrons to be captured by each indi- 
vidual nucleus. After freeze-out, there are two important 
facts to be considered: 1) nuclei compete for capturing 
the few neutrons available, 2) the neutron-capture and 
beta-decay rates are comparable. These two facts are 
common to the hot and cold r-process and are important 
for determining the final abundances shown in Fig. [5] by 
the green lines. Here we focus on the difference among 
mass models and general features will be described in the 
next section. 

Figure [5] shows the abundances obtained using the four 
mass models for the hot (left column) and cold (right col- 
umn) r-process. The freeze-out abundances (black lines) 
are characterized, especially in the hot r-process, by the 
presence of strong fluctuations that have almost disap- 
peared in the final abundances (green lines), as expected 
from solar r-process abundances. These fluctuations are 
due to the fact that (n, 7) and (7,^) reactions favor nu- 
clei with an even neutron number. Consequently, as Z 
increases in moving from one isotopic chain to the next, 
N increases by at least two units (except at the magic 
numbers where it stays constant). Therefore, some mass 
numbers are not present in the r-process path as shown 
by the dots in Fig. [6] 

The freeze-out abundances can be understood looking 
at the two neutron separation energy in Fig. [6l Two 
kind of features in S^n leave a fingerprint on the abun- 
dances: 1) The abrupt drop in S^n at the magic numbers 
N = 82 and N — 126, leads to accumulation of matter 



at these neutron numbers and to the formation of peaks 
in the abundance distribution. 2) In regions where S2 n is 
constant or presents a saddle point behaviour, an equi- 
librium can not be achieved between neutron captures 
and photodissociation. In the hot r-process, this leads 
to troughs in the abundances around A ~ 110,140 for 
the FRDM mass model and around A ~ 185 for FRDM, 
ETFSI-Q, and HFB-17 mass models. These features of 
the two-neutron separation energies have less impact for 
the cold r-process because photodissociation reactions 
are suppressed due to the low temperatures. 

We use the average neutron separation energy and the 
neutron-to-seed ratio (see Fig. [7|) to discuss the evolution 
of matter during the r-process for different mass models 
in the hot r-process. The average neutron separation 
energy, (S n ), shows two minima and one maximum for all 
mass models. However, (S n ) significantly differs during 
the early evolution in the calculation based on ETFSI- 
Q nuclear masses. For t « 30 ms, when matter flow 
approaches N = 82 shell closure, (S n ) is smaller and the 
minimum is broader. Due to the quenching of the N = 82 
shell gap introduced in the ETFSI-Q mass model [6l|, the 
abrupt drop in S2n at N — 82 dissapears for Z < 43 (see 
Fig. [6j. Consequently, with this mass model the r-process 
proceeds through nuclei with smaller neutron separation 
energy before reaching the N — 82 shell closure. This 
occurs at Z = 43 for ETFSI-Q, while for the other mass 
models it is reached already for Z = 40. The width 
of the first minimum is related to the sum of half-lives 
of the nuclei on the r-process path with TV < 82. This 
is substantially larger in the calculations with ETFSI-Q 
mass model. There is thus a delay in the time required to 
overcome the N = 82 shell closure and a slow down of the 
speed at which neutrons are captured. Therefore, the r- 
process freeze-out occurs at later times for ETFSI-Q than 
for HFB-17 and Duflo-Zuker. The situation is different 
for FRDM. Here the neutron separation energy drops 
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FIG. 4: (Color online) Neutron density, neutron-to-seed ratio, 
and average neutron separation energy evolution (Eq. ((3j ) for 
the same trajectories shown in Fig. Q] 



peaks at 138 Sn in the hot r-process calculation and at 
134 Cd and 140 Sn in the cold r-process (see upper panels 
of Fig. [5]). These nuclei represent a barrier to the flow 
of neutron captures to more neutron-rich isotopes and 
consequently the matter has to wait for their beta-decay 
before heavier nuclei can be reached. This effect results in 
a broader second minima for the FRDM average neutron 
separation energy in Fig. [7] and in a longer duration for 
the r-process. 

In the phase after freeze-out, the few available neu- 
trons are not equally captured in all regions. This leads 
to different position of troughs and peaks depending on 
the mass model used. The most visible feature in the 
final abundances is the trough in the ETFSI-Q abun- 
dances before the third peak. This is also present in the 
freeze-out abundances based on FRDM and HFB-17, but 
not on Duflo-Zuker. During the decay to stability the 
trough is filled when using the FRDM and HFB-17 mass 
models while for the ETFSI-Q becomes even larger. The 
two neutron separation energies in Fig. [5] show that for 
FRDM there is a drop of S 2n for N = 122 followed by a 
rise before the N = 126 magic shell. This leads to the for- 
mation of the trough at A ~ 184. For ETFSI-Q the two 
neutron separation energies are almost constant for nuclei 
in the region A = 180-190 and Z ~ 60. This produces 
two troughs in the freeze-out abundances at A = 180 
and A = 187. The situation is more complicated for 
HFB-17 (as the neutron separation energies show larger 
fluctuations) with the net result of troughs at A = 180 
and A = 184 in the freeze-out abundances. During the 
decay to stability, in the calculations based on FRDM, 
HFB-17, and Duflo-Zuker neutron captures move matter 
from the region before the trough to higher mass numbers 
and the trough is partially filled. In contrast, ETFSI-Q 
presents higher S2n (and thus higher (71,7) rates) in the 
region just before N — 126. This leads to a shift of mat- 
ter from the trough towards the peak that produces an 
enhancement of the first one. 

Notice that Duflo-Zuker abundances present stronger 
odd-even effects than the abundances calculated with the 
other models as shown in bottom panels of Fig. [5] How- 
ever, this effect is not due to the mass model itself but 
to the computed neutron-capture rates, which here are 
based on the simple approximation suggested in Ref. [67| . 
The importance of neutron-capture rates on the final 
abundance will be discussed in detail in Sects. IHICI and 
MICH 



C. Decay to stability 



abruptly just before N = 90 and even becomes negative 
for nuclei like 133 Pd, 134 Ag, and 137 Cd. This region is 
reached when the r-process breaks out of the N — 82 shell 
closure. As the beta-decay half-lives of these nuclei are 
relatively long (around 100 ms for the rates used in the 
present calculation |70j ]) matter accumulates producing 



As we have shown in the previous section, there are still 
reactions occurring after freeze out that contribute to the 
redistribution of matter and to the production of the final 
abundances. Classical r-process studies (see for example 
Ref. [HI, but still amply used for r-process chronome- 
ters studies [74], [Zlj]) neglect the neutron captures after 
freeze out and consider that beta-delayed neutron emis- 
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1. The hot r-process path is shown at freeze-out (Y n /Y Bee d = 1) by dots. 



sion is the only mechanism to redistribute and smooth 
the abundances after freeze-out. Dynamical r-process 
calculations (26j have shown that neutron captures dur- 
ing freeze-out can reduce odd-even effects but also shift 
the peaks [3lj and produce the small rare-earth peak 
around A rj 160 |3(J. However, in some of these studies 
(e.g., HEIHj]) neutron captures were suppressed once the 
Y n /Y see d ratio was small and only beta-decays were con- 
sidered. Here, we show that even when Y n /Y sec d « 10~ 5 , 
neutron captures are still key for determining the final 
abundances. For these conditions the neutron density 
is around 10 18 cm~ 3 , leading to a typical time between 
neutron captures of 100 ms, that is comparable with the 
beta-decay half-lives. 

The processes producing the smoothing of the r- 
process abundances act mainly between freeze out 
(Y n /Y scc( \ = 1) and the moment when Ti n ,^) = T /3- The 
upper panels in Fig. [8] show the abundances at freeze- 
out (black lines) and when the average neutron cap- 
ture (Eq. (p?aj) ) and beta-decay (Eq. Pc|0 timescales are 
identical (green lines) for hot (left column) and cold 



(right column) r-process. The abundances at freeze- 
out present large fluctuations (particularly in the hot re- 
process) which have almost completely disappeared at 
the later time as indicated by the green line. There is 
still some redistribution of matter at even later times 
that leads to the final abundances shown in Fig. [5] and 
to the formation of the rare-earth peak around A rs 160. 

The competition between beta decay (with delayed- 
neutron emission) and neutron capture when matter de- 
cays to stability is important to understand the differ- 
ences in the final abundances between the hot and cold 
r-process. Let us consider a nucleus with charge Z and 
mass number A. Its abundance can change by beta- 
decay, neutron capture, and photodissociation. The com- 
petition among these processes can be quantified by the 
beta decay flux: 



Fp{Z,A) = X p {Z,A)Y{Z,A) ; 



(4) 



10 




FIG. 7: (Color online) Evolution of the average neutron separation energy and neutron-to-seed ratio for the hot (left column) 
and cold (right column) r-process and based on the four mass models discussed in the text. 



and by the net neutron capture flux: 

F n (Z,A) = Y(Z,A)N n (av} z , A -X J {Z,A+l)Y(Z,A+l). 

(5) 

In order to visualize these quantities, it is convenient 
to define the fluxes for an isotopic chain, Fp(Z) = 
Ea F p(Z, A ) and F n (Z) = "£ A F n {Z,A), and the 
fluxes for an isobaric chain, Fp(A) = Fp(Z, A) and 
F n (A) = ^ z F n {Z 1 A). As discussed in the Appendix, 
we expect that for the hot r-process a beta-flow equi- 
librium is achieved (23[. This implies that the Fp(Z) 
reaches a constant value independent of Z. In the cold r- 
process one expects that both Fp{Z) and F n (A) become 
constant. This is confirmed in Fig. [5] that shows the 
net neutron capture and beta-decay fluxes versus mass 
number (middle panels) and versus proton number (bot- 
tom panels). Notice, that in the s-process [76| F n (A) 
is also constant and there are strong odd-even effects in 
the abundances. This is due to the large odd-even effects 
present in the neutron capture rates and the fact that 
beta-decay can be assumed instantaneous compared to 
s-process timescales. In contrast, in our case these strong 



odd-even effects are not present because beta-decay rates 
become similar to neutron-capture rates as matter decays 
to stability. 

The fluxes F n (A) and Fp{A) present several features 
when Tp = T( n 7 ) that can explain how matter is redis- 
tributed. In the regions where beta-decay dominates over 
neutron captures, nuclei will beta decay without sub- 
stantially changing the mass number, e.g., see Fig. [5] 
for A > 195. On the other hand, nuclei in regions 
where neutron capture dominates over beta-decay will 
predominantly capture neutrons and consequently the 
abundances will shift to higher mass numbers, as shown 
in Fig. |5] for A = 182-195. The formation of the rare- 
earth peak (not yet present in the abundances shown in 
upper panels of Fig. [8]) is also due to neutron capture as 
matter decays to stability. In the region A — 162-168 the 
beta-decay and neutron-capture fluxes are very similar, 
while in the region A < 162 the latter dominates (Fig. [5]). 
This produces a net movement of matter from nuclei with 
A < 162 to nuclei with A w 162 that will result in the 
formation of the rare-earth peak. The formation of this 
peak has to wait until the moment when the beta decay 
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fluxes become larger than the neutron capture fluxes in 
that region. The reason why the neutron capture fluxes 
can become locally smaller for these nuclei (when com- 
pared with slightly heavier or lighter nuclei) is the pres- 
ence of a deformed sub-shell closure around A w 162 (see 



Fig. IH]) that results in a sudden drop of neutron-capture 
rates. 

In the hot r-process, the fluxes (Fig. [5J right column) 
present more fluctuations due to photodissociation reac- 
tions. These reactions can result in negative net neutron 
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capture fluxes in regions where the photodissociation is 
still important. This leads to an extra supply of neu- 
trons that can be important at later times. Negative net 
neutron capture fluxes appear in the hot r-process for 
A < 145 and Z < 53 once Tp — T(„ )7 ). They are not 
shown in Fig. |8]as we use a logarithmic scale. 



1. The role of neutron capture 

In order to explore the impact of neutron capture we 
compare two different sets of rates based both on the 
FRDM mass model. The first set corresponds to statis- 
tical model calculations of Ref. [64| and it was used in 
previous sections. The second set is com put ed with the 
analytic approximation suggested in Ref. |67| . These two 
sets of neutron-capture rates are compared in Fig. [9] for 
three different isotopic chains in regions relevant for the 
r-process. The Ru isotopes are populated in the region 
N < 82, Xe isotopes for N ~ 100, and Er isotopes for 
N < 126. The largest differences between both sets of 
neutron-capture rates occur for nuclei just before neutron 
shell closures. For these nuclei the level density around 
neutron separation energy becomes rather low and con- 
sequently the rates are very sensitive to the treatment of 
parity [77| and to the dipole strength distribution [78$ . 
In addition, the statistical model may not be applica- 
ble for some of these nuclei at r-process temperatures 
(see Fi g. 7 of Ref. [7^) and direct capture should be in- 
cluded [H, Ho| . 




Neutron Number N 

FIG. 9: (Color online) Neutron-capture rates as computed in 
Ref. [6^ ] (solid lines and filled symbols) compared to rates 
based on the analytical approximation of Ref. (dashed 
lines and empty symbols). The rates are shown for 3 dif- 
ferent isotopic chains in the region of relevance for r-process 
nucleosy nt hesis . 

The abundances based on the two sets of neutron- 
capture rates are shown in Fig. [10] for the hot and cold r- 
process. At freeze-out both sets of neutron-capture rates 
leads to very similar abundances. In the hot r-process, 
the abundances are independent of the neutron-capture 



rates since the evolution proceeds under (n, 7) T± (7, n) 
equilibrium (see Appendix). The cold r-process is char- 
acterized by the competition between neutron capture 
and beta decay, however only small changes are present 
in the freeze-out abundances when the neutron-capture 
rates are varied. Notice that the position and the height 
of the peaks are the same. In contrast, the final abun- 
dances exhibit significant differences: The third r-process 
peak is more shifted towards higher mass number and the 
abundances between peaks are less smooth for the calcu- 
lations based on the approximate neutron-capture rates. 

The competition of the nuclei to capture the few neu- 
trons available can be quantified by the net probability 
of a nucleus for neutron capture, that we define as 

using net neutron capture flux, F n (Z,A), introduced in 
Eq. ([5])- The denominator in this expression represents 
the change of neutron abundance due to neutron cap- 
tures and photodissociations and it is always positive 
since the neutron abundance continuously decreases dur- 
ing an r-process calculation. The numerator is positive 
for the majority of nuclei, although can become nega- 
tive for some of them. The change in abundances after 
freeze-out is thus more pronounced in regions with larger 

Pn,y{Z, A). 

Figure [IT] shows the neutron-capture probability for 
different regions: around the second r-process peak (A < 
140), between peaks (140 < A < 185), and around the 
third peak (A > 185) for the hot r-process. The results 
are qualitatively the same for the cold r-process. As ex- 
pected from the agreement between the freeze-out abun- 
dances, the initial evolution of the probabilities is rather 
similar and follows the build up of increasingly heavier 
nuclei during the r-process. At early times most of the 
captures takes place in the region around and below the 
second r-process peak. Later as nuclei in the region be- 
tween peaks are produced, the neutrons are mainly cap- 
tured in this region. Just before freeze-out (t » 1 s), 
there is an increase in the capture probability around 
the third peak. The evolution after freeze-out strongly 
depends on the neutron captures. In the calculation with 
NON-SMOKER rates (Fig. [TIJ left panel), the neutron- 
capture probability is similar for the regions around the 
third peak and between peaks, with the latter dominating 
slightly. This is in contrast to the neutron-capture prob- 
ability based on the approximate rates (Fig. [TTJ right 
panel) which is clearly higher in the region of the third 
peak than in the region between peaks. Therefore, neu- 
trons are mainly captured around the third peak. This 
leads to the shift of the peak towards even higher A than 
with the NON-SMOKER rates and leaves pronounced 
odd-even features between peaks. Consequently, when 
the approximate rates are used, beta-delayed neutron 
emission is the only mechanism to smooth the abun- 
dances between A = 130 and A = 195 after freeze out. 
While using the NON-SMOKER rates, both (neutron 
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FIG. 10: (Color online) Abundances at freeze-out (upper panels) and after decay to stability (bottom panels) for the hot (left 
panels) and cold (right panels) r-process. The abundances are obtained using two different sets of neutron-capture rates both 
based in the same nuclear mass model (FRDM [6(|). The first set of neutron-capture rates, labeled NON-SMOKER (black 
lines), corresponds to the calculations of Ref. [13]. The second set, labeled as "Approximation" (green lines), has been obtained 
using the analytical approximation derived in Ref. [63]. 



capture and beta-delayed neutron emission) contribute 
to produce a smother r-process distribution. 

Our results demonstrate that neutron-capture rates are 
key for the determination of the final r-process abun- 
dances, specially around the third peak and between 
peaks where a robust pattern is found in old metal-poor 
halo stars and solar system (see e.g., Ref. [zl]). Simi- 
lar conclusions were reached in Ref. [81[ where individ- 
ual neutron-capture rates were changed by an arbitrary 
factor. However, our results illustrate more clearly the 
non-local character of the competition for the few avail- 
able neutrons which can produce global changes in the 
r-process abundances. In addition to the variation of the 
third-peak position, we find also substantial changes in 
the abundances around A ~ 150, even if both regions are 
not directly connected by any reaction. Similar global 
changes were obtained in Refs. [2!| [32|, [82| where the 
sensitivity of the r-process abundances to changes of the 
neutron-capture rates around A = 130 was studied. 



2. Beta-delayed neutron emission 

After we have shown the importance of beta decay 
when matter decays to stability, it is worth to analyze the 
effect of beta-delayed neutron emission. In the classical 
r-process, where no neutron captures are considered after 
freeze out, beta-delayed neutron emission is the only way 
to redistribute matter and get smooth final abundances. 
In Fig. [TH we explore the effect of beta-delayed neutron 
emission in the hot (left column) and cold (right column) 
r-process. The black lines correspond to the standard 
network calculations with beta-delayed neutron emission 
included and the green lines to calculations where beta- 
decay takes place without emitting neutrons, i.e. A is 
conserved. The effect of beta-delayed neutron emission 
and its subsequent capture was also investigated for dif- 
ferent long-time dynamical evolutions in [29l . |38| . 

After freeze-out, more neutrons are present when beta- 
delayed neutron emission is considered and thus the 
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FIG. 11: (Color online) Probability of neutron capture (Eq. ([6])) vs. time for different mass number intervals based on NON- 
SMOKER rates (left panel) and on the approximated rates (right panel). 



neutron-to-seed ratio is higher as shown in the bottom 
panels of Fig. Q21 However, in our hot r-process such 
differences are too small to have an impact on the fi- 
nal abundances. In this case, photodissociation produces 
also neutrons in addition to beta decay. Our results seem 
to contradict classical r-process calculations, where the 
redistribution of matter is due only to beta-delayed neu- 
tron emission. One should notice that our freeze-out path 
is in agreement with classical r-process calculations (see 
e.g. [23]). Moreover, if we neglect all neutron captures 
after freeze-out (as it is done in the classical r-process) 
and only consider beta-delayed neutron emission, we ob- 
tain the abundances shown by the red line in Fig. [T2] 
This calculation reproduces the position and width of 
the third peak but present larger oscillations in the fi- 
nal abundances. Similar results were also found in the 
classical r-process calculations of Ref. [23[ demonstrating 
that beta-delayed neutron emission cannot completely re- 
move the fluctuations in the freeze-out abundances. Once 
neutron captures are considered the abundance distribu- 
tion becomes smoother like in the solar system and the 
rare-earth peak forms (black line in the upper panels of 
Fig. [T2)) . However, the third peak becomes narrower and 
shifts to mass number values larger than A ~ 195. 

In our cold r-process, the beta-delayed neutron emis- 
sion has more impact on the final abundances (Fig. IT2|) . 
Since photodissociation is negligible, the r-process path 
can move farther away from stability reaching nuclei with 
higher probability of emitting neutrons after beta-decay. 
This leads to significant differences in the neutron-to-seed 
ratio shown in the bottom panel of Fig. [12"J The freeze- 
out is more instantaneous without beta-delayed neutron 
emission and this has two main effects: the third peak 
is less shifted and the rare earth peak is not produced. 
Therefore, we can conclude that beta-delayed neutron 
emission is very important to supply the neutrons that 
through several captures will determine the final abun- 



dances. Moreover, the freeze-out cannot be totally in- 
stantaneous because neutron capture are required to form 
the rare earth peak at late times. 



3. Non-instantaneous freeze-out 

Finally, we want to discuss a general feature present in 
the evolution of the neutron-to-seed ratio. After an ini- 
tial slow decrease, the neutron-to-seed reaches a phase of 
fast decline once its value becomes around one. How- 
ever, this fast decline, that will correspond to an in- 
stantaneous freeze-out as assumed in classical r-process 
calculations, is always interrupted and the neutron-to- 
seed ratio follows a more moderate decrease afterward. 
This is a generic feature found in all dynamical calcu- 
lations (see e.g. Refs. [lg, |25|, |83[) and indicates that 
the freeze-out effects discussed here will always be im- 
portant. We found that the sudden change in the evo- 
lution of the neutron-to-seed ratio, mathematically cor- 
responding to the appearance of an inflexion point, oc- 
curs when the average neutron-capture rate (Eq. (|2a[) ) 
and the average beta-decay rate (Eq. ([2c]) ) are identi- 
cal, i.e. T(„ 7 ) = Tp. Interestingly, this happens even in 
the calculations where beta-delayed neutron emission is 
artificially switched off, see Fig. [12"J This phase of mod- 
erate decline of the neutron-to-seed ratio is the so-called 
s-process phase of the r-process [84[ in which beta-decay 
dominates over neutron-capture. During this phase, the 
rare-earth peak is formed in our calculations based on 
the FRDM mass model. However, in these calculations 
the third r-process peak is shifted to mass number values 
larger than A ~ 195. 
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FIG. 12: (Color online) Effect of the beta-delayed neutron emission in the hot (left column) and cold (right column) r-process, 
using neutron captures and beta decays based on the FRDM mass model. The black lines are for the reference case calculated 
with the standard nuclear input: neutrons are emitted with given probability (P n ) after beta decay. The green lines are for the 
case P n — 0, i.e. A is conserved during beta decay. For the red dash-dotted line the neutron captures and the photodissociation 
reactions are suppressed after freeze-out. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the impact of the long-time dynami- 
cal evolution and nuclear physics input on the r-process 
abundances. Our calculations are based on hydrodynam- 
ical trajectories from core-collapse supernova simulations 
of Ref. |9j with the entropy increased by a factor two (i.e., 
density decreased by a factor two) in order to produce the 
third r-process peak. We have chosen two different evo- 
lutions to cover the two possible physical conditions at 
which the r-process occurs in high entropy ejecta. These 
evolutions are identical during the seed formation phase 
and differ only after the temperature becomes < 3 GK 
and the nucleosynthesis flow is dominated by neutron 
captures, i.e. the r-process phase. This guarantees that 
changes in the resulting abundances are due to the mod- 
ification of the long-term evolution and/or the nuclear 
physics input and not to changes in the initial condi- 
tions. The long-time evolution is varied assuming that 



the reverse shock is at different temperatures, which is 
justified based on two-dimensional simulations [46j. The 
two typical long-time evolutions are: 

• hot r-process that occurs when the reverse shock 
is at high enough temperatures (T > 0.5 GK) to 
reach an (n, 7) (7, n) equilibrium that lasts until 
neutrons are exhausted. 

• cold r-process (with the reverse shock at low tem- 
peratures) that takes place under a competition be- 
tween neutron capture and beta decay. 

The main difference between these evolutions is that in 
the cold r-process the photodissociation is negligible [39j . 
Therefore, the r-process path can move farther away from 
stability reaching nuclei with shorter beta-decay half- 
lives and leading to a faster evolution and an earlier freeze 
out than in the hot r-process. 
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We can distinguish two phases during the r-process 
that are characterized by different nuclear physics pro- 
cesses. Before freeze-out, Y n /Y sec <i > 1, the most relevant 
nuclear physics input depends on the type of r-process 
taking place, hot or cold r-process: 

• In the hot r-process the most relevant input are the 
nuclear masses as they determine the r-process path 
via the neutron separation energy. A very good ap- 
proximation to the freeze-out abundances of the hot 
r-process is obtained assuming that matter achieves 
steady beta- flow, i.e. \r{Z)Y{Z) = constant (see 
Eq. (|A.7jl ). It is well known [85, 86[ that in this case 
the peaks in the r-process abundance distribution 
at A = 130 and A — 195 are associated with long 
beta-decay half-lives at the magic neutron numbers 
N = 82 and N = 126 where the r-process path is 
closer to the stability. 

The abundances at freeze-out in the hot r-process 
are characterized by the presence of large troughs 
that occur in regions where the two neutron sepa- 
ration energy is almost constant or presents a sad- 
dle point. This typically occurs before (after) a 
magic neutron number where a transition from de- 
formed (spherical) to spherical (deformed) shapes 
takes place. As experimental data for two neutron 
separation energies do not show this behavior [87l | , 
this may indicate a drawback in some of the mass 
models used. However, the solar abundances sug- 
gest a small trough in the region before the third 
peak, A = 180-190, that could be related with a 
transition from deformed to spherical nuclei. 

• In the cold r-process the most relevant inputs are 
beta-decay and neutron-capture rates. The nu- 
clei at the r-process path are those with similar 
neutron-capture and beta-decay rates for a given 
neutron density. Our results are the first to show 
that the abundances at freeze-out achieve steady 
flow for both beta decays and neutron captures, 
i.e. Xp(Z)Y(Z) and (<jv)aY(A) are constant (see 
Eq. (IA.8I) 1 ). This suggests that the cold r-process is 
more robust than the hot r-process, as the abun- 
dances fulfil additional constraints. 

The sensitivity of the mass model has been investi- 
gated by consistently using neutron separation energies 
and neutron capture rates based on the mass models: 
FRDM, ETFSI-Q, HFB-17, and Duflo-Zuker. Our re- 
sults show peculiarities coming from each mass model. 
In ETFSI-Q the quenching of the N = 82 shell closure 
leads to a slow down of the evolution and to a later freeze- 
out. Moreover, the large values of S2n before N = 126 in 
this mass model make the trough in the freeze-out abun- 
dances for A rs 185 bigger due to neutron captures when 
matter decays to stability. Results based on FRDM are 
clearly affected by the anomalous behaviour of S2 n be- 
fore N = 90, which produces the accumulation of matter 



and thus the formation of peaks around A w 135 even in 
the cold r-process. 

In order to study the evolution after freeze-out we have 
used the fluxes for neutron captures and beta decays. 
They help us to explain the final features in the abun- 
dances, such as the exact position of the r-process peaks 
and the formation of the rare earth peak. These are our 
most significance outcomes: 

• The abundances at freeze-out can be approximated 
assuming a steady beta-flow (hot r-process) or a 
steady flow of beta decays and neutron captures 
(cold r-process) for a given neutron density and 
temperature. However, the final abundances are 
determined by the evolution after freeze-out. In 
all cases considered, the final abundances are sub- 
stantially different and smother than the freeze-out 
abundances. Most of the smoothing takes place 
just after freeze-out when the timescale for neu- 
tron captures is still shorter than the one for beta- 
decays, T(„ j7 ) < Tp, see Eq. (J2|). Hence, neu- 
tron captures play a dominant role in producing 
a smooth distribution. Nevertheless, the evolution 
once T7„ 7 ) > Tp is also important. During this 
phase the decrease of the neutron-to-seed ratio is 
rather moderate and determined by the timescale 
at which matter beta decays to stability, even in 
calculations where beta-delayed neutron emission 
is artificially suppressed. 

• The impact of the neutron capture rates has been 
investigated by comparing results based on the 
same mass model (FRDM) but different sets of 
neutron-capture rates: A set is based in sta- 
tistical model calculations with the code NON- 
SMOKER [64] and the other in the analytical ap- 
proximation derived in Ref. [67]. We find that the 
abundances at freeze-out for hot and cold r-process 
are rather similar for both sets of neutron capture 
rates. However, after freeze-out we find that most 
of the neutrons are captured in the region between 
r-process peaks when using the NON-SMOKER 
rates. In contrast, with the approximated rates 
neutrons are captured more probably in the region 
around the third peak. The end result is a larger 
shift of the third peak and a less smooth abundance 
distribution between peaks with the approximated 
rates than with the NON-SMOKER rates. This 
emphasizes the important role of neutron captures 
after freeze-out. 

• The small rare-earth peak, observed around A ~ 
160 in the solar r-process distribution, must neces- 
sarily be formed after the freeze-out, since it is not 
present in any of the freeze-out abundances. Fur- 
thermore, it is also not present in the abundances 
when T(„ i7 ) w Tp. We find, see also Refs. (30l. l3lj|. 
that the rare-earth peak forms by neutron captures 
when matter decays to stability. 
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• The main role of beta-delayed neutron emission is 
the supply of neutrons. In the hot r-process we 
find no difference in the abundances calculated with 
and without beta-delayed neutron emission. Since 
temperature is high, photodissociation prevents the 
path from reaching regions far from stability where 
the probability of emitting neutrons after beta de- 
cay is large. Furthermore, photodissociation reac- 
tions are also a source of neutrons. In contrast, in 
the cold r-process the neutron-to-seed ratio reaches 
significantly smaller values when beta-delayed neu- 
tron emission is suppressed. This leads to a reduc- 
tion in the shift of the third peak after freeze out 
but also inhibits the formation of the rare earth 
peak. This confirms the argument given above that 
the peak forms by neutron captures. 

Our results clearly rise the importance of future ex- 
periments to measure nuclear masses, neutron capture 
rates and beta-decay half-lives for nuclei far from sta- 
bility. This will provide not only direct input for net- 
work calculations, but also important constraints for the 
theoretical nuclear models. We have shown that the r- 
process abundances are very sensitive to the set of neu- 



tron capture rates as they determine the regions in which 
neutrons are capture predominantly. More experimental 
effort is necessary for an improved determination of neu- 
tron capture cross sections. Since these experiments are 
difficult, sensitivity studies to determine the most rele- 
vant neutron capture rates will be necessary. Our results 
show a strong interplay between the late-time evolution 
of the ejected matter and the nuclear physics input. This 
could constrain the astrophysical conditions once future 
radioactive experimental facilities deliver high quality ex- 
perimental data for r-process nuclei. 



Appendix: Formalism 

During the r-process, and assuming that charged parti- 
cle reactions, fission and alpha decays can be neglected, 
the evolution of the abundances is mainly determined 
by neutron capture, photodissociation, and beta-decays. 
This results in the following differential equation that de- 
termines the change of the abundance of a nucleus with 
charge Z and mass number A: 



J 



dY(Z,A) 
dt 



pN A (av) z , A -iY n Y(Z, A - 1) + A 7 (Z, A + l)Y(Z, A + l) 
j 

+ Xpjn{Z - 1, A + j)Y(Z -1,A + j) 



pN A (av} z , A Y n + X 7 (Z,A)+J2^n(Z,A) ] Y{Z,A) 

3=0 



(A.l) 



r 



where Y n is the neutron abundance, {av(Z, A)} is the neutrons (up to a maximum of J). The photodissociation 

thermal averaged neutron-capture rate, and A 7 (Z, A) the rate is related to the neutron capture rate by detailed 

photodissociation rate for a nucleus Z, while \pj n (Z, A) balance: 
is the f}~ decay rate of A Z with emission of j delayed 



J 



\{Z,A + 1) = (av) 



Z.A 



,.kT 



2irh 2 



3/2 



2G(Z,A) ( A 



G(Z,A + 1) \ A + l 



r 



3/2 



exp 



S n (Z,A+l) 



kT 



(A.2) 



where G is the particition function and S n = m n + 
M(Z, A — 1) — M(Z, A) is the neutron separation energy 
with m n the neutron mass and M(Z,A) the mass of the 
nucleus. 

If the assumption is made that the neutron abundance 
varies slowly enough, it can be assumed that the neu- 



tron density, N n = Y n pN A , is constant over a timestep. 
In this case the network can be divided into separate 
pieces for each isotopic chain and solve then sequentially, 
beginning with the lowest Z 56|. However, this approxi- 
mation becomes numerically unstable when the neutron 
abundance becomes small, Y n < 1CP 5 . Consequently, it 
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is better to include in the set of differential equations the 
one determining the change of the neutron abundance: 



dt 



Y,pNA(<Jv) z , A Y n Y{Z,A) 



Z,A 



J2\(Z,A)Y(Z,A) 



+ Y.i X ^n{Z,A) Y(Z,A) (A.3) 

Z,A \j=l J 

The system of differential equations defined by 
Eq. (jA.ip and (|A.3|) allows for several approximations 
that are valid in different physical regimes. A commonly 



used assumption in classical r-process calculations is the 
(ri,7) ^ (7,") equilibrium. This approximation is valid 
whenever the neutron density (N n > 10 20 cm" 3 ) and 
temperature (T > 1 GK) [8g] are large enough to war- 
rant that both the rate of neutron capture (N n (av) and 
the photodissociation rate (A 7 ) are much larger than the 
beta decay rate (Xp) for all the nuclei participating in 
the network. Under this conditions the evolution of the 
system is mainly determined by the beta decay rates 
as the abundances along an isotopic chain are inmedi- 
atly adjusted to an equilibrium between neutron cap- 
tures and photodissociations, i.e. N n Y(Z, A){av) z,A — 
X 7 ( Z, A + 1)Y(Z, A + 1). Combining these result with 
Eq. ([A~2]) one obtains that the abundances in an isotopic 
chain are given by the simple relation: 



Y(Z,A + 1) 
Y(Z,A) 



2nh 2 
m u kT 



3/2 



A + l \ 3/2 G(Z,A + 1 ) 
A J 2G(Z, A) 



r 



exp 



S n {Z.A + l) 



kT 



(A.4) 



For each isotopic chain, the above equation defines a nu- 
cleus that has the maximum abundance and which is 
normally known as waiting point nucleus as the flow of 
neutron captures "waits" for this nucleus to beta-decay. 
The set of waiting point nuclei constitutes the r-process 
path. The maximum of the abundance distribution can 
be determined setting the left-hand side of Eq. (|A.4|) to 
1, which results in a value of S n that is the same for all 
isotopic chains for a given neutron density and tempera- 
ture: 



S°(MeV) 



T 9 
5.04 



34.075 - log N n + - log T g , (A.5) 



where Tg is the temperature in units of 10 9 K and N n is 
the neutron density in cm -3 . Equation (|A.5|) implies that 
the r-process proceeds along lines of constant neutron 
separation energies towards heavy nuclei. For typical r- 
process conditions this corresponds to 5° ~ 2-3 MeV. 
Due to pairing, the most abundance isotopes have always 
an even neutron number. For this reason, it may be more 
appropriate to characterize the most abundance isotope 
in an isotopic chain as having a two-neutron separation 
energy S^n — 25," [8^ |. The two- neutron separation en- 
ergy is not a continuous function of neutron number but 
shows large jumps particularly close to magic neutron 
numbers. For this reason r-process nuclei near to magic 
numbers have neutron separation energies much larger 
than the typical 2-3 MeV and the r-process path moves 
closer to the stability (see figure [6]). 

If (n, 7) (7, n) equilibrium is valid it is sufficient to 
consider the time evolution of the total abundance of an 
isotopic chain Y(Z) = J^a y {Z, A) as the abundances 



of different isotopes are fully determined by Eq. (|A.4|) . 
From Eq. (|A.1|) we can determine the time evolution of 
Y(Z) obtaining: 



dY(Z) 
dt 



Xp{Z -\)Y{Z -\)- Xp{Z)Y{Z) (A.6) 



where Xp(Z) = £ A A^Z, A)Y(Z, A)/Y(£). In this 
case the r-process evolution is independent of the 
neutron-capture rates, only beta-decays are necessary for 
Eq. (1X6] and masses via S n in Eq. (|A.4[) . If the r-process 
proceeds in (n, 7) (7, n) equilibrium and its duration is 
larger than the beta decay lifetimes of the nuclei present, 
Eq. (|A.6I) tries to reach an equilibrium denoted as steady 
/3-fl.ow [23J that satisfies for each Z value: 



Xp(Z-l)Y{Z-l)=X p (Z)Y{Z) 



(A.7) 



In this case the peaks at A = 130 and 195 in the solar r- 
process distribution can be attributed to the long /3-decay 
lifetimes of the waiting point nuclei with N — 82 and 126, 
where the r-process path gets closer to the stability (see 
figure ■ This is the case for the equilibrium calcula- 
tions discussed in the text before freeze-out of neutron 
captures. 

The r-process can also operate under such a low tem- 
peratures that the photodissociation rates in Eq. (|A.1|) 
can be neglected. Under this conditions the r-process 
operates under a competition of neutron captures and 
beta decays. If one neglects beta-delayed neutron emis- 
sion, Eq. (|A.1[) can be reduced to two independent equa- 
tions that govern the evolution of the total abundance 
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along an isotopic chain, Y(Z) and along an isobaric chain, 
Y(A)=J2 Z Y(Z,A): 



dY{Z) 



= Xp(Z-l)Y(Z-l)-Xp{Z)Y(Z) (A. 



dY ( A ) »T / I 



-iY(A — 1) — N n ((jv)AY(A) (A.8b) 

were {ctv)a = 'Y!, z {ov)z,aY(Z,A)/Y{A). If the r-process 
duration is longer than the beta-decay and neutron cap- 
ture lifetimes, Eq. (|A.8[) reaches an equilibrium that we 
will denote as steady flow that satisfies for each Z and 
A: 



X P {Z -1)Y(Z -1) = X (Z)Y{Z) (A.9a) 



(av) A . 1 Y(A - 1) = (av) A Y{A) (A.9b) 

In addition, as the r-process occurs under a competition 
of beta-decays and neutron captures one obtains that 



N n (crv)AY(A) w Xp(Z)Y(Z). As the abundances along 
an isotopic and isobaric chain are dominated by a sin- 
gle nucleus this condition determines also the nuclei that 
participate in the r-process, i.e. the r-process path. Sim- 
ilarly to what happens in the equilibrium case, the peaks 
in the abundance distribution correspond to long beta 
decay lifetimes. However, in this case the peaks are in 
addition associated with long neutron capture lifetimes. 
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